Anomalous Dynamic Arrest in a Mixture of Big and Small Particles 
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We present molecular dynamics simulations on the slow dynamics of a mixture of big and small 
soft-spheres with a large size disparity. Dynamics are investigated in a broad range of temperature 
and mixture composition. As a consequence of large size disparity, big and small particles exhibit 
very different relaxation times. As previously reported for simple models of short-ranged attractive 
colloids and polymer blends, several anomalous dynamic features are observed: i) sublinear behavior 
for mean squared displacements, ii) concave-to-convex crossover for density-density correlators, by 
varying temperature or wavevector, iii) logarithmic decay for specific wavevectors of density-density 
correlators. These anomalous features are observed over time intervals extending up to four decades, 
and strongly resemble predictions of the Mode Coupling Theory (MCT) for state points close to 
higher-order MCT transitions, which originate from the competition between different mechanisms 
for dynamic arrest. For the big particles we suggest competition between soft-sphere repulsion 
and depletion effects induced by neighboring small particles. For the small particles we suggest 
competition between bulk-like dynamics and confinement, respectively induced by neighboring small 
particles and by the slow matrix of big particles. By increasing the size disparity, a new relaxation 
scenario arises for the small particles. Self-correlators decay to zero at temperatures where density- 
density correlations are frozen. The behavior of the latters resembles features characteristic of 
type-A MCT transitions, defined by a zero value of the critical non-ergodicity parameter. 

PACS numbers: 82.70.Dd, 64.70Pf, 83.10.Rs 



I. INTRODUCTION 

Rheological properties of soft matter based systems 
can be manipulated by a proper addition of com- 
ponents of different mobilities. Some examples are 
colloid-polymer mixtures and polymer blends. These 
type of binary mixtures exhibit unusual relaxation fea- 
tures which challenge standard pictures for structural 
dynamic arrest in glass-forming liquids and colloidal 
systems. Differently from the usual two-step increase 
and decay for, respectively, mean squared displace- 
ments and dynamic correlators [ij, |j, yi Oj ^] ? the lat- 
ters do not exhibit a defined plateau at intermediate 
times between the microscopic and diffusive regimes 
0, 0, lal- This result suggests a softer character for 
the collective caging mechanism — i.e., the temporary 
trapping of each particle by its neighbors. Dynamic 
correlators show a highly stretched decay, and for se- 
lected values of the control parameters the decay is 
logarithmic in time. By varying wavevectors or con- 
trol parameters the decay shows a striking concave- 
to-convex crossover [^ Q , leading to violation of stan- 
dard scaling laws for complex dynamics as the time- 
temperature superposition principle. 

These anomalous relaxation features have been re- 
cently derived within the Mode Coupling Theory 
(MCT) for simple models of hard-sphere colloids with 
short-ranged attractions [9j. The MCT is a first- 
principle theory of the glass transition which makes 
predictions for averaged dynamic quantities as mean 
squared displacements, diffusivities or incoherent and 
coherent dynamic correlators [ij, lillO]. The only in- 
put in MCT equations is the structural information 
contained in total and partial static structure factors, 
which enter the memory function accounting for the 
fluctuating forces. MCT determines dynamic quan- 



tities from the knowledge of static correlations. In- 
deed MCT relates small variations in static correla- 
tions with large variations in dynamics. 

Solution of the MCT equations for the mentioned 
models of short-ranged attractive colloids [9|, [ill, [l^ 
I13I Il4j has determined the existence of so-called 
higher-order transitions, showing properties rather 
different from the standard fold MCT transition as- 
sociated to the liquid-glass transition. The mentioned 
anomalous relaxation features have been related to 
the presence of nearby higher-order transitions. MCT 
predictions for these type of models have been con- 
firmed by molecular dynamics simulations |a, 13, Il5| 
and experiments [l^ |l3, ll8| . 

The models of hard sphere-colloids with short- 
ranged attractions investigated in Refs. [a, 0, 13 are 
used as one-component effective models for the col- 
loidal particles in colloid-polymer mixtures. Addi- 
tion of small polymers (or other small particles) to 
dense solutions of big colloidal particles yields an ef- 
fective attraction between the colloidal particles in or- 
der to maximize entropy |l9l |. This effect is known 
as the depletion mechanism. As a consequence, in 
a certain range of density, temperature and mixture 
composition, competition occurs between two differ- 
ent mechanisms for dynamic arrest of the colloidal 
particles: hard sphere-repulsion characteristic of col- 
loidal systems, and formation of reversible bonds, in- 
duced by the depletion mechanism |20|, |2l| . In the ef- 
fective one-component system the higher-order MCT 
scenario arises as a consequence of these two compet- 
ing mechanisms of very different localization lengths 
9, 13]. When heating up or cooling down the system, 
dynamic arrest is exclusively driven by, respectively, 
hard-sphere repulsion and reversible bond formation, 
and relaxation features of standard liquid-glass tran- 



sitions are recovered 

Very recently, we have carried out simulations on 
a simple bead-spring model for polymer blends 8]. 
The introduction of a signifficant monomer size dis- 
parity yields very different relaxation times for both 
components in the blend, the component of small 
monomers being the fast one. The fast component 
exhibits anomalous relaxation features, very differ- 
ent from standard results observed for homopolymcrs, 
and strongly resembling predictions of the higher- 
order MCT scenario for short-ranged attractive col- 
loids. Fully atomistic simulations on a real poly- 
mer blend are consistent with the anomalous features 
reported for the bead-spring model |22|. We have 
pointed out the hypothesis of an underlying higher- 
order MCT scenario for the dynamics of the fast com- 
ponent, which might arise from the competiton be- 
tween bulk-like caging and confincmcmt as different 
mechanisms for dynamic arrest for the fast compo- 
nent. Bulk- like caging is induced by the particles of 
the fast component, and confinement is induced by the 
matrix formed by the chains of the slow component. 
Due to chain connectivity, the former mechanism is 
present even for high dilution of the fast component, 
extending the anomalous relaxation scenario over a 
broad range of blend compositions IM] . 

At present no MCT theoretical calculations are 
available for models ofpolymer blends similar to that 
investigated in Ref. j^]. A related system for which 
theoretical MCT works are available is a binary mix- 
ture of big and small hard spheres |23|. However, the 
size disparity used in such works, and also in computa- 
tional investigations in the MCT framework 24, f2M > 
is not sufficient to provide a large separation in the 
time scales for big and small particles, i.e., to induce 
confinement effects for the latters. It is not also suffi- 
cient to induce signifficant depletion effects for the big 
particles. As a consequence, the anomalous relaxation 
features reported for polymer blends and short-ranged 
attractive colloids are not predicted for moderate size 
disparity, and indeed are not observed in the corre- 
sponding simulations. Instead, a description in terms 
of standard M CT p redictions is possible for moderate 
size disparities iHElEsf. 

Very recently, MCT theoretical calculations have 
been reported by Krakoviack |26ll27l | for a binary mix- 
ture of mobile and static hard sphereSja system where 
confinement effects are present ,23, 1221 • For this sys- 
tem, the extreme cases of the Lorentz gas (a single 
mobile particle in a disordered medium of static ob- 
stacles) and the liquid of hard spheres are obtained for 
high dilution of respectively mobile and static parti- 
cles. At a given composition of the mixture a higher- 
order MCT transition has been derived. This result 
supports the hypothesis of a similar MCT scenario in 
polymer blends originating from competition between 
bulk- like caging and confinement. 

The observed analogies between dynamics of col- 
loidal particles in colloid-polymer mixtures, and of 
the fast component in polymer blends suggest that the 
higher-order MCT scenario might be a general feature 



of systems showing slow dynamics with several com- 
peting mechanisms for dynamic arrest. In this article 
we provide new evidence in favour of this hypothesis 
by carrying out molecular dynamics simulations on a 
mixture of big and small soft-spheres of very different 
sizes. Results presented here for large size disparity 
(ratio S = 2.5) complement previous investigations on 
slow dynamics in binary mixtures with small disparity. 
We observe anomalous relaxation features similar to 
those recently reported for models of short-ranged at- 
tractive colloids and polymer blends. We investigate 
a wide range of temperatures and mixture composi- 
tions. By tuning the composition, these features are 
displayed both by the big and the small particles. 

For the case of the big particles, we assign such 
anomalous features to competition between soft- 
sphere repulsion and the depletion mechanism induced 
by the small particles. For the small particles, we sug- 
gest competition between bulk-like dynamics induced 
by the neighboring small particles, and confinement 
induced by the matrix of slow big particles. Sim- 
ilarly to the fast component in polymer blends |3], 
and despite of the absence of chain connectivity, small 
particles exhibit apparent anomalous relaxation over 
a broad range of compositions extending up to high 
dilution. In the latter case, these effects are clearly 
manifested by a signifficant subset of small particles 
forming small clusters. 

We have also performed simulations at a fixed com- 
position for very large size disparity {d = 8.0). A new 
relaxation scenario arises for the small particles, show- 
ing features characteristic of nearby MCT transitions 
of the so called type-A. Such transitions are defined by 
a zero value of the long-time limit (non-ergodicity pa- 
rameter) of density-density correlators, different from 
the finite value defining the usual type-B transitions. 
This feature provides a connection with results in Ref. 
l2Q| for a mixture of mobile and static particles, which 
report a dynamic phase diagram displaying an A- and 
a B-line merging at a higher-order point. Hence, we 
suggest that results at moderate disparity 5 = 2.5 
for the system here investigated might originate from 
the existence of a nearby B-line (providing finite val- 
ues for the non-ergodicity parameters) ending at a 
nearby higher-order point (providing anomalous re- 
laxation features). 

The article is organized as follows. In Section II 
we introduce the investigated model and give compu- 
tational details. In Section III we present simulation 
results for static structure factors. We also present 
dynamic quantities displaying unusual relaxation fea- 
tures. In Section IV the framework of the MCT is 
used in an operational way to describe simulation re- 
sults. In Section V we discuss the possible origin of 
the observed anomalous dynamic features. We also 
propose a picture for the different relaxation scenar- 
ios observed by increasing the size disparity of the 
particles. Conclusions are given in Section VI. 

II. MODEL AND SIMULATION DETAILS 

We have simulated a mixture of big (labelled as A 



and B) and small (C and D) particles of equal mass 
TO = 1, interacting through a soft-sphere potential 
plus a quadratic term: 
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where e = 1 and a, (3 e {A, B, C, D}. The inter- 
action is zero beyond a cutoff distance ccjap, with 
c = 1.15. The addition of the quadratic term to the 
soft-sphere interaction, with the values Cq = 7c~^^ 
and C2 = Gc"^"*, guarantees continuity of potential 
and forces at the cutoff distance. The diameters of the 
soft-sphere potential for the different types of interac- 
tion are: ctdd = 1, crcc = I-Ictdd, o-bb = 2.5tTDD, 
CTAA — I-Icbb, and CTq^ = {(Jaa + o'/3/3)/2 for the case 
a^f3. 

The potential (^ is purely repulsive. It does not 
show local minima within the interaction range r < 
caap- Hence, slow dynamics in the present model 
arises as a consequence of steric effects. MCT the- 
oretical works are usually carried out on systems of 
hard objects, while simulations in similar systems with 
continuous interactions are usually preferred for com- 
putational simplicity. In the present system, the tail 
of the interaction potential is progressively probed by 
decreasing temperature, which plays the role of in- 
creasing packing in a system of hard spheres. Hence, 
simulations presented here should be useful for a qual- 
itative test on the success or failure of future MCT 
theoretical works on mixtures of hard spheres with 
large size disparity. 

The composition of the mixture is defined as the 
fraction of small particles: Xsmaix — {Nq + N^) / {Na + 
Nb + Nc + N-£,), with Na denoting the number of 
particles of the species a. As shown below, the in- 
troduction of a large size disparity between the sets 
{A,B} and {C,D} yields very different time scales for 
both sets. We impose the constraints A'a = Nb and 
Nq ~ Nu- These constraints, together with the small 
selected ratios crcc/cDD — caa/cbb = 1.1 guarantee 
that only very small dynamic differences are induced 
between particles within a same set ({A,B} or {C,D}), 
and at the same time, avoid crystallization for the in- 
vestigated compositions Xsmaii- Crystallization would 
occur for very asymmetric mixtures if onl y on e type 
of big and small particles were introduced |2fl| • 

In the following, temperature T, distance, wavevec- 
tor q, and time t, will be given respectively in units 
of e/kB, Ctdd, '^bd^ ^^'^ crDD("^/e)^'^. The packing 
fraction, 6, is defined as: 
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+Ncalc + A^d^d] (2) 

with L the side of the simulation box. Simulations 
have been carried out at a constant packing frac- 
tion (f) = 0.53. This value is comparable to those 
used in simulations of slow relaxation in simple liq- 
uids. For comparison, the original Lennard- Jones bi- 
nary mixture investigated by Kob and Andersen [3j 



has (j) = 0.59 with the definition of packing frac- 
tion given above. We investigate the T-dependence 
of the dynamics for mixture compositions a;smaii = 
0.1, 0.3, 0.6 and 0.8. The number of big and small 
particles for each composition are respectively {Na + 
Nb:Nc + Ad)= (5400:600), (2100:900), (1000:1500) 
and (800:3200). The system is prepared by placing the 
particles randomly in the simulation box, with a con- 
straint that avoids core overlapping. Periodic bound- 
ary conditions are implemented. Equations of motion 
are integrated by using the velocity Verlet scheme [31| , 
with a time step ranging from 2 x 10""* to 5 x 10~^, 
for respectively the highest and the lowest investigated 
temperature. A link-cell method |3]J is used for saving 
computational time in the determination of particles 
within the cutoff distance of a given one. 

At each state point, the system is thermalized at the 
requested temperature by periodic velocity rescaling. 
After reaching equilibrium, energy and pressure show 
no drift. Likewise, mean squared displacements and 
dynamic correlators show no aging, i.e., no time shift 
when being evaluated for progressively longer time ori- 
gins. Once the system is equilibrated, a microcanoni- 
cal run is performed for production of configurations, 
from which static structure factors, mean squared dis- 
placements, and dynamic correlators are computed. 
For each state point, the latter quantities are aver- 
aged over typically 20-40 independent samples. 

III. RESULTS 

a. Static structure factors 

We compute normalized partial static structure 
factors, Sapiq) = {Pa{q, 0)p/3(-g, 0))/y/NaNf3, with 
Pail, t) = Sj exp[iq • raj(i)], the sum extending over 
all the particles of the species a G {A,B,C,D}. Fig. 
n] shows, at a fixed temperature T = 0.50 and differ- 
ent mixture compositions, results for A-A, D-D, and 
A-D pairs. Data for other big-big, small-small, and 
big-small pairs display only small quantitative differ- 
ences. 

For high concentration of the big particles, Xsmaii = 
0.1, S'aa(<z) shows a sharp first peak at g = 2.65a^^, 
corresponding to a typical distance of 2.37crDD be- 
tween neighboring A-particles. This distance is 
smaller than the soft-sphere diameter ctaa = 2.75crDD- 
This feature is possible due to the interpenetrable 
character of soft spheres. By increasing the concen- 
tration of small particles, the matrix of big particles 
progressively becomes more disordered. As a conse- 
quence, the height of the first peak of Saa{q) decreases 
considerably. The first minimum follows the opposite 
trend and becomes less pronounced. The small shift of 
maxima and minima to higher g- values by increasing 
a^smaii reflects a stronger packing of the big particles. 
For large concentration of small particles, Xsmaii = 0.8, 
'S'aa(Q') shows a small peak at low-g. We assign this 
peak to the presence of inhomogeneities or "voids" 
in the matrix of big particles, which are filled by the 
small particles (see also Section V, Fig. lilt;'). 

At high dilution of the small particles, Xsmaii = 0.1, 
Sub{,<i) shows a nearly structureless profile, close to 
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FIG. 1: (color online) Partial static structure factors at 
T = 0.50 and different mixture compositions, as computed 
from simulation data. Wavevectors are given in units of 



the flat behavior expected for a gas. Only a weak os- 
cillation is observed at \ow-q. This feature reflects the 
existence of a small fraction of clusters of neighbor- 
ing small particles (see also Section V). Increasing the 
concentration of small particles yields a broad peak 
at g w G.ScTpQ, corresponding to a typical distance of 
about (TDD between neighboring small particles. The 
peak grows up and narrows by increasing Xsmaih as 
a signature of progressive ordering of the small parti- 
cles. Still, one finds a rather broad peak at the highest 
value of a;smaii = 0.8. Increasing the concentration of 
small particles also produces a low-g peak of increas- 
ing intensity in Suuiq)- Such a peak originates from 
the inhomogeneities in the structure formed by the 
small particles (see also Section V, Figs. ITTb andlTTl:). 
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FIG. 2: (color online) For the investigated compositions, 
temperature dependence of the diffusivities for the dif- 
ferent species of the mixture. For low temperatures at 
a^smaii = 0.1, data for the big particles are absent, since 
they do not reach the diffusive regime in the time window 
of the simulation (see text and Fig. |3^) . 



We do not observe changes, at any of the investigated 
temperatures, in the intensity of the low-q peak within 
the tinre window of the simulation. Hence, phase sep- 
aration is discarded for the results here presented. 

b. Diflfusivities and mean squared 
displacements 

Fig. [3 shows the T-dependence of the diffusivity, 
D, for all the species in the mixture at all the investi- 
gated compositions. For each species a € {A,B,C,D} 
the diffusivity is calculated from the long time limit 
of {[Ara{t)f)/6t, with {[Arait)f} the corresponding 
mean squared displacement at tinre t. As shown in 
Fig. [3 the introduction of a signifficant size disparity 
yields very different time scales for the sets of big and 
small particles. Small differences are instead obtained 
between the diffusivities of both species within a sanre 
set. Only for the lowest temperature at concentration 
a^smaii = 0.1 there are signifficant differences (about a 
factor 4 in diffusivity) between C- and D-particles. In 
the following, simulation results will only be shown for 
the big A-particles and the small D-ones. The quali- 
tive behavior of respectively B- and C-particles is the 
same, displaying only small quantitative differences 
with the formers. 

Fig. Oshows the T-dependence of the mean squared 
displacement of the A-particles, ((ArA)^}, for three 
different compositions Xsmaii — 0.1, 0.6 and 0.8. As 
usually observed in the proximity of liquid-glass tran- 
sitions 0) 0) 13 ) a bending occurs after the initial bal- 
listic (oc i^) regime. A plateau arises at low tem- 
peratures. This effect corresponds to the well-known 
caging regime — i.e., the temporary trapping of each 
particle in the cage formed by the neighboring ones. 
At long times, the diffusive regime (oc t) is reached 
for values ((ArA)^} ^ "'aA' i-^' when the A-particles 
have moved, on average, a distance of the order of 
their size. 
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FIG. 3: (color online) Symbols: mean squared displace- 
ment of A-particles for different compositions of the mix- 
ture. Straight lines indicate linear or subhnear behavior. 



By looking in more detail at the data for T — 0.75 
and T = 1.0, at Xsmaii > 0.6 (Figs. Eb and Et), 
an unusual appproximate sublinear behavior, ex i", 
is observed over two time decades after the ballistic 
regime. The exponent a < 1 decreases by decreas- 
ing temperature. For ccsmaii — 0.8 we also note, at all 
temperatures, a bump at the interval 0.2 ^ t < 0.4. 

Fig. ^ shows, at the same compositions, results 
for D-particles. Remarkable differences with mean 
squared displacements of A-particles are observed. 
Differently from A-particles (and from the standard 
behavior in the proximity of liquid-glass transitions) 
D-particles reach the diffusive regime oc t for displace- 
ments much larger than their size: from ((ArD)^) 
IOct^d at high T to ((Aro)^) 



100ct|)d at low T. 



FIG. 4: (color online) As Fig. I^for D-particles. 



This result is observed for all the compositions. As in 
the case of the A-particles, for Xsmaii > 0.6 an unusual 
approximate sublinear regime is observed at interme- 
diate times, with an exponent decreasing by decreas- 
ing temperature. This sublinear regime sets on for 

Results reported in this subsection evidence the ex- 
istence of unusual relaxation features in the slow dy- 
namics of mixtures of big and small particles with 
sufficiently large size disparity. Next we evaluate the 
effects of size disparity in partial density-density cor- 
relators. 

c. Density-density correlators 

We compute partial density-density correlators, 
Faf3{q,t) = {pa{q,t)pi3{-q,0))/{pa{q,0)pi3{-q,0)). 
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FIG. 5: (color online) Symbols: T-dependence of 
FAA{q,t), at specific q- values, for several compositions. 
The straight line in panel (c) indicates logarithmic decay. 

Fig. shows the T-dependence of FAA{q,t) at spe- 
cific values of the wavevector q and different compo- 
sitions. For the case ccsmaii = 0.1 (Fig. |5Jl), A- A 
correlations display the standard behavior observed 
for Hquid-glass transitions 0, S S Hi B- After the 
initial transient regime, Faa(9,0 shows a first decay 
to a plateau. By decreasing temperature, the plateau 
extends over longer time intervals. At long times, a 
second decay occurs from the plateau to zero. This 
second decay corresponds to the a-process of the glass 
transition and is usually well described by a stretched 
exponential function. 

Figs. IHId and[Sl: show results for the case Xsmaii = 
0.8, for wavelengths {2n/q) of respectively 2.2ctdd and 
I.Ictdd- While apparently standard behavior is ob- 



FIG. 6: (color online) As Fig. Elfor FDD(g,t). 



tained in the former case, unusual features are ob- 
served for wavelengths probing the size of the small 
particles (Fig. Et). First, Faa(9, i) does not exhibit 
a defined plateau. Moreover, the shape of the long- 
time decay shows a concave-to-convex crossover by 
decreasing temperature. At an intermediate temper- 
ature r = 0.50 the decay is logarithmic over two time 
decades. As in the case of the mean squared displace- 
ment, a bump is observed at the interval 0.2 < t ^ 0.4. 

Unusual relaxation features are also displayed by 
FDD{q,t). Fig. El shows results for different composi- 
tions at wavelengths probing the size of the big par- 
ticles. No defined plateaux are observed for ajsmaii = 
0.1 and 0.6. Logarithmic relaxation is observed at 
intermediate temperatures for time intervals extend- 
ing up to three decades. Data for Xsmaii = 0.6 



exhibit a concave- to-convex crossover. Due to the 
mentioned partial crystalhzation of big particles for 
a^smaii = 0.1 [301 , data are absent in Fig. |S^ for the 
range 0.50 < T < 1.0. Apparently, a concave-to- 
convex crossover is also present for this composition, 
though data are not conclusive. 

Results for i^DoC?, t) at Xsmaii = 0.8 display a quali- 
tatively different behavior. At the lowest investigated 
temperature T = 0.30 a bump is observed at the in- 
terval 3 < t < 6, followed by a decay until t sa 60, 
where logarithmic relaxation sets on and extends over 
three time decades. 

IV. MCT ANALYSIS 

a) Big particles 

Many of the anomalous relaxation features reported 
in Section 111 — sublinear behavior for mean squared 
displacements, and logarithmic decay and concave- 
to-convex crossover for density-density correlators — 
strongly resemble those reported for hard-sphere col- 
loids with short-ranged attractions 6, 7], and for the 
fast component in polymer blends with components 
of very different mobilities [Sl, |22| . As mentioned in 
the Introduction, these anomalous features have been 
strictly derived in the framework of the MCT for sim- 
ple models of short-ranged attractive colloids J9| . Ac- 
cording to MCT, this anomalous relaxation scenario 
arises from an underlying higher-order transition. Mo- 
tivated by this fact, we discuss the present results by 
using MCT in an operational way. 

In its ideal version, which neglects activated hop- 
ping events, MCT predicts a sharp transition from 
an ergodic liquid to a non-ergodic arrested state at 
a given value Xc of the relevant controlparameter x 
(in practice density or temperature) [ij, ^ [lOj . When 
crossing the transition point x = Xc from the er- 
godic to the arrested state, the long-time limit of the 
density-density correlator for wavevector q, F{q,t), 
jumps from zero to a non-zero value, denoted as the 
critical non-ergodicity parameter, /^. Moving beyond 
the transition point into the non-ergodic state yields a 
progressive increase of the non-ergodicity parameter, 
fq > fq- In the MCT formalism, the standard liquid- 
glass transition is of the fold type, also denoted as A2 
in Q> llO yj| ' In the standard case the jump in F{q, t) 
is discontinuous, i.e., the critical non-ergodicity pa- 
rameter /^ takes a finite value. MCT transitions with 
/^ > are also denoted as type-B transitions. For 
ergodic states close to the transition point, the ini- 
tial part of the a-process — i.e., the von Schweidler 
regime — is approximated by a power law expansion 



Fiq,t)^r^-hqit/r)'' + h(^\tM 



26 



(3) 



.(2) 



with < b < 1. The prefactors hq and h\ ' only 
depend on the wavevector q. The characteristic time 
T only depends on the separation parameter \x — Xc\ 
and is divergent at the transition point. Another im- 
portant prediction of the MCT for state points close 



to fold transitions is the so-called second universality 
or x-time superposition principle (with x the corre- 
sponding control parameter). According to this pre- 
diction, the final decay of F{q,t) (i.e., the final part 
of the a-process) is unvariant under scaling by the a- 
relaxation time Ta- Hence, for two state points x — xi 
and X = X2 close to the transition point the final decay 
of F{q, t) fulfills the relation 



(4) 



tion_pomt 



Fx=xi q, 



Ta{xi) 



= Fa:=x2 q, 



Taix2) 



The a-relaxation time Ta is a time scale probing the 
a-process. In practice, it can be obtained from fitting 
the decay from the plateau to a stretched exponential 
Aq exp[—(t/Ta)^], with Aq the plateau height and < 
/3 < 1. It can also be defined as the time where F{q, t) 
decays to some small value, e.g. 0.3 (tq.s), provided 
it is well below the plateau. 

Another prediction of the MCT for points close to 
an A2-transition is the power-law decrease of the dif- 
fusivity to zero, D ex |a; — xd'^ . The exponent 7 is 
given by 0,1 Q: 



7 



1 
2a 



1 
26' 



(5) 



with < a < 0.4. Hence 7 > 1.75. The critical expo- 
nents a, b and 7 are related to the so-called exponent 
parameter 1/2 < A < 1 through the relation: 



A 



T^l + b) _ T^{l-a) 
r(l + 26) " r(l-2a)' 



(6) 



with r the Gamma function [ij, |3, llOf . 

By numerically solving the MCT equations for dy- 
namic correlators, transition points are determined as 
those where the respective long time limit exhibits 
a jump from a zero to a non-zero value. From the 
knowledge of the total and partial static correlations 
at the transition point, all the critical exponents and 
the coefficients in Eqs. (|21 El El a-re univoquely de- 
termined 0, y, llOl- Solving the MCT equations is 
a difficult task, which in general is only possible for 
rather simplified models of real systems. Hence, in- 
stead of solving the equations, critical exponents and 
prefactors are often obtained as fit parameters from 
simulation or experimental data. Consistency of the 
data analysis requires that the so-obtained critical ex- 
ponents fulfill Eqs. (|3J) and ©. 

Fig. 13 shows an analyisis of data for A-particles 
at composition Xsmaii = 0.6 in the framework of the 
MCT for A2-transitions. Fig. [7|i shows a fit of the 
diffusivity to a power law Da oc (T — Tc)'^. The fit 
provides the values Tc — 0.180 and 7 = 3.80. From 
Eqs. © and ® we obtain A = 0.887, a = 0.220, and 
b = 0.326. Fig. {7}p shows, for FaaC?, t) at a low T, fits 
of the decay from the plateau to Eq. ||3J). By forcing 
the von Schweidler exponent to the mentioned value 
b = 0.326 a consistent description is achieved. Fig. 
Eh shows a succesful test of the second universality 
(time-temperature superposition principle), by using 
To. 3 for the a-relaxation time. For the composition 
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FIG. 7: (color online) MCT analysis of data for A- 
particles, for composition a^smaii = 0.6. Panel (a): T- 
dependence of the diffusivity. Symbols are simulation 
data. The line is a fit to a MCT power law Da oc (T-Tc)^ 
(see text). Panel (b): Symbols: g-dependence of FAA(g, i) 
at a low temperature. Curves are fits to Eq. I^. Panel 
(c): T-dependence of FAA{q,t) for a specific g- value. In 
order to test the second universality, times are divided by 
the a-relaxation time ro.3(r) (see text). 



a^smaii = 0.3 similar results are obtained (not shown). 
From an analogous analysis we obtain T^ = 0.344, 
A = 0.880, 7 = 3.68, a = 0.226, and b = 0.340 for this 
latter composition. 

For the composition Xsmaii = 0.1, a test of MCT 
for A-particles is not possible due to the absence of 
data in the range 0.50 <T < 1.0. Though a two-step 
decay is observed for T < 0.50 (see Fig. |3Ji), a fit 
to Eq. © is not possible in that temperature range, 



FIG. 8; (color online) Panel (a): Symbols: g-dependence 
of FAA{q,t) for a composition a;smaii ~ 0.8 at a low tem- 
perature. Curves are fits to Eq. ||7J. Panels (b) and (c): 
Values of, respectively, /^ and Hq, as obtained from fits 
to Eq. of FAA{q,t) at two state points SI and S2. The 
characteristic times r for SI and S2 are respectively 9 and 
3. 



since the plateau height is not constant, but clearly 
increases by decreasing temperature. According to 
predictions of ideal MCT, these temperatures would 
correspond to non-ergodic states below the transition 
point (T < Tc). Within this interpretation, the fact 
that a decay from the plateau occurs at long times, 
would be a signature of the so-called hopping events, 
which restore ergodicity. 

For the composition Xsmaii = 0.8, data of FAAi,q,t) 
for wavelengths 27r/g probing the size of the small par- 
ticles do not show a defined plateau (Fig. \^), and a 



fit to Eq. (O is not possible. Moreover, the observed 
logarithmic decay and concave-to-convex crossover do 
not fit to expectations for A2-transitions. They in- 
stead resemble the features reported for short-ranged 
attractive colloids |6|,|7| and polymer blends ^, which 
for the former case have been derived in the framework 
of MCT for state points close to higher-order transi- 
tions. 

A higher-order MCT transition, A„+i, is character- 
ized by an exponent parameter A = 1 , and can emerge 
as the result from the interplay between n control 
parameters {x"} = {xi,X2---Xn)- Higher-order MCT 
transitions were initially derived for schematic mod- 
els 123, UM I but only recently have been obtained for 
short-ranged attractive colloids [^ |lj [l^, llM UM ^^ 
a first realization in real systems. It can be demon- 
strated that an expansion in a power-law series as Eq. 
is not convergent for A = 1 or for values of A 
very close to unity. On the contrary, a logarithmic 
expansion is rapidly convergent. Hence, close to a 
higher-order MCT transition, or more generally to an 
A2-transition with A ^, 1, F(ci,t) is approximated at 
intermediate times as [sl, l22l bJi 1^^ '■ 

Fiq, t) « r^ - H, ln(t/r) + i/f ) InHt/r), (7) 

(2) 

where the prefactors Hq and Hq depend on q and on 
the distance of the state point {x"} to the transition 
point {x"}. These prefactors exhibit two important 
properties Q: ''^) Hq factorizes as Hq — C{{s."-})Hq 
where Hq only depends on q, and the ^-independent 
term Cdx"}) depends on the state point. Hence, val- 
ues of Hq for different state points close to the tran- 
sition point must be proportional, ii) H\ does not 
follow scaling behavior. It is decomposed as |9J: 



il(2) = ^({x"}) -f B({x"})if, 



«' 



(8) 



where Kq only depends on q, and the g-indcpendent 
terms yl({x"}) and i?({x"}) depend on the state 
point. The terms Hq and Kq are univoquely deter- 
mined by static structure factors at the higher-order 
transition point. As in the case of A2-transitions, coef- 
ficients in Eq. Q are often obtained as fit parameters 
from simulation or experimental data. 

Decomposition of H\ according to Eq. JHJ has an 
important consequence. There are hypersurfaces in 
the control parameter space, {x"} = {x"}((7), where 
Hq changes its sign, being zero along the hypersur- 
face. This property leads, for a given value of q, to 
a concave-to-convex crossover in F{q, t) when cross- 
ing the hypersurface by varying control parameters 
(as temperature or density) '9J. Analogously, for a 
given state point, varying the value of q also leads to a 
concave-to-convex crossover in F{q, t). Since Hq ' — 
for state points at the hypersurface, according to Eq. 
O F{q,t) will exhibit a logarithmic decay for such 
state points. Moving between different state points 
changes the value of q for which pure logarithmic de- 
cay occurs. The concave-to-convex crossover is one of 
the main signatures of the higher-order MCT scenario 



and differentiates it from other theoretical frameworks 

i 

As shown in Figs. \^ and|SjL, a concave-to-convex 

crossover is present, both by varying temperature and 
wavevector, for FAA{<l,t) at composition a;smaU = 0.8. 
Fig. Hi shows fits to Eq. for T == 0.50 M- 
Figs. IHt andlHt show respectively the obtained val- 
ues for fq and Hq, for temperatures T = 0.40 and 
T = 0.50. The fact that a common f^ is found for 
both states, together with the observed scaling behav- 
ior of Hq, would be consistent with the existence of a 
nearby MCT higher-order transition (A = 1) or an A2- 
transition with A ^ 1. Differently from the usual be- 
havior for liquid-glass transitions, fq shows no strong 
modulation but a nearly monotonous decay. Though 
this result qualitatively resembles the observed q- 
dependence at A3- and A4- transitions for models 
of short-ranged attractive colloids [a, y, |l3, llJ, I^TJ ]. 
such a comparison can be misleading. The observed 
g-dependence is not necessarily related to a hypothet- 
ical higher-order MCT transition. It is indeed also 
found at moderate size disparit y, fo r which standard 
relaxation features are observed |23, |2a| . 

It is also noteworthy that the observed sublinear 
behavior in the mean squared displacement at inter- 
mediate times is another feature characteristic of the 
higher-order MCT scenario. In leading-order [2j, the 
mean squared displacement is given by ([Ar(i)]^) = 
r1{t/TY, with Tc the localization length and z = 
^MSDC'({x"})/r^. The coefficient /imsd is determined 
by static structure factors at the transition point. 
Since the prefactor C({x"}) decreases as the tran- 
sition point is approached, the exponent z also de- 
creases, as observed in Fig. 13:. Deviations from 
pure power-law behavior result from corrections to the 
leading term. Such corrections are minimal for certain 
values of the control parameters, yielding pure power- 
law behavior |2|. 

b) Small particles 

Data for D-particles cannot be reproduced by the 
standard A2-scenario for any of the investigated com- 
positions, due to the absence of defined plateaux and 
to the presence of logarithmic relaxation in density- 
density correlators. Data for Xsmaii = 0.6 exhibit fea- 
tures resembling predictions of the higher-order MCT 
scenario. Fig. El shows results from such an analy- 
sis, analogous to that in Fig. |S1 for A-particles, in 
terms of the latter framework. The concave-to-convex 
crossover in Fdd(9, t) obtained by decreasing temper- 
ature at constant q (Fig. Eb), is also observed by vary- 
ing q at constant temperature (Fig. E^l). A good de- 
scription of the decay is achieved by Eq. O . Consis- 
tently, fits for different state points provide the same 
values for fq (Fig. EJs), and scaling behavior is ob- 
tained for the corresponding values of Hq (Fig. Et). 
The observed sublinear behavior for the mean squared 
displacement, and the decrease of the exponent by 
decreasing temperature (Fig. ^p) is another feature 
shared with the higher-order MCT scenario. 

We have performed an analogous analysis for 



10 



nO.8 

Q 

Q 
Ph 0.6^ 



0.4 

0.2 





rHk!Ks^"'s3: X ' ,, = 0.6' T = 03 

^^i^'^t-^ small 


. q=1.38 
q = 2.09 




- q = 2.45 
> q = 2.66 
. q = 2.93 


X, ■ * N; -^X ^^ V "v 


► q = 3.22 


xX^\XXS^v\ 


+ q = 4.01 


\X* XXtvxX 


» q = 4.46 




. q = 5.20 
. q = 6.27 




(a) ^^""""^-^SS^S^Sl " 


Mil 1 1 1 1 1 



10 



10 



lo' 10^ 



3 4 

10 t 10 



1- 

0.8- 

IT 
I 

0.6- 

0.4- 

0.2- 

0, 







0.1 



0.05- 



S3: X ,, = 0.6 T = 0.3 

small 

S4: X ,, = 0.6 T = 0.4 

small 



5. 



O fq(S3) 

n fq(S4) 



(b) 



Q 9. 



e @ 



- 9 

(c) 




S3: 
S4: 




X 
X 

9 


sma„ = 0-6 T = 0.3 
sma„ = 0-6 T = 0.4 




O 1.26H^(S3) 
n H^(S4) 




s 


= «B 






10"' lo" lo' 10^ 



0.8- 



0.6 



0.4 



0.2 



0.1- 
0.08- 
0.06 
0.04 
0.02 



3 4 

10 t 10 



- n 
o 


S5: 
S6: 


small 
small 


= 0.1 T = 0.3 
= 0.1 T = 0.4 


8 








O fq(S5) 




- 


a 

D 






n fq(S6) 


_ 






- 


(b) 

1 


1 


1 


1 1 1 1 



"0 1 







: 9 

' (c) 


8 


Q 


S5: 
S6: 


small 
small 


= 0.1 T = 
= 0.1 T = 


0.3 _ 
0.4 




QH (S5) 
D H (86) 


- 


^ 


'segg 


) - 







FIG. 9: (color online) As Fig. Elfor FoD(q,t), x^^m = 0.6, 
at the state points S3 and S4. The characteristic times r 
for S3 and S4 are respectively 30 and 7.2. Data in panel 
(a) for the wavevector q = 2.45(Tq]3 exhibit logarithmic 
decay over four time decades. 



FBDiq,t) at composition Xsmaii — 0.1 (see Fig. [TU|) . 
despite there are some differences with features char- 
acteristic of the higher-order MCT scenario. Differ- 
ently from the case Xsmaii = 0.6, no apparent subhnear 
behavior is present in the mean squared displacement 
(Fig. ^). Due to the absence of data in the range 
0.5 < T < 1.0, it is difficult to unambiguously identify 
a concave-to-convex crossover in -Fdd {q, t) by varying 
temperature (Fig. Et). Fig. ITUk shows results by 
varying q at constant T = 0.30. The trend exhibited 
by the data, which show a extremily stretched decay 
at large wavevectors, suggests that such a crossover 
might be present for higher values of q. However, we 



FIG. 10: (color online) As Fig. IHlfor Fdd ((?,*), 3;sniaii = 
0.1, at the state points S5 and S6. The characteristic times 
r for S5 and S6 are respectively 360 and 50. 



cannot confirm this point, which remains to be under- 
stood. We have not detected the crossover at least up 
to q Ki ^CTpQ. Beyond that q- value the amplitude of 
the decay is rather small and it is difficult to solve the 
shape of the curve within statistical noise. 

Still, Eq. (TJ provides a good description of the 
decay (Fig. ilOb). with a common f^ for different state 
points (Fig. llOb). and Hq obeying scaling behavior 
(Fig. llOb). Surprisingly, the scaling factor is unity 
within error bars, despite dynamics at state points 
T — 0.30 and T — 0.40 are signifficantly different (see 
Figs. Hl and Ell). 

Due to the complex form of the decay of i^DD(<Z, t) 
for composition Xsmaii = 0.8 (Fig. Et), a fit to Eqs. 
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© or lO is not possible. Hence, at this composi- 
tion there are no apparent analogies for small particles 
with known MCT scenarios. 

V. DISCUSSION 

Results presented in Sections III and IV exhibit 
strong dynamic analogies with MCT predictions for 
state points close to higher-order transitions (A = 1), 
or more generally, for A2-transitions with exponent 
parameter A ^ 1. These analogies must not be un- 
derstood as a proof of an underlying MCT scenario 
for the mixture of big and small particles here in- 
vestigated. An unambiguous answer to this ques- 
tion could only be provided by solving the corre- 
sponding MCT equations. As mentioned in the In- 
troduction, theoretical and computational works on 
the framework of the MCT on mixtures of big and 
small hard spheres have not reported anomalous re- 
laxation features |2J, |2J, |2^. However these works 
have explored size disparities smaller than the value 
S = (Taa/ccc = cbb/cdd = 2.5 used in this work. 

Fig. ^3 displays typical slabs of the simulation box 
for mixture compositions Xsmaii — 0.1, 0.6, and 0.8. 
The slab thickness is 5(Tdd ■ On the basis of particular 
features displayed by the configuration of small and 
big particles at the different compositions, next we 
discuss the observed anomalous relaxation features in 
terms of competition between different arrest mecha- 



a) Symmetric mixtures 

For composition Xsmaii = 0.6, big particles are dis- 
tributed over the simulation box in a rather homo- 
geneous way, with a weak tendency to form clusters. 
Small particles fill connected cavities in the slow ma- 
trix of big particles, which acts as a confining medium 
for the formers. As a consequence, competition be- 
tween bulk-like dynamics and confinement occurs for 
arrest of the small particles. Pure logaritmic decay in 
FBD{q,t) is observed for q ~ 2.5a^l) (Fig. |5^), i.e., 
for wavelengths ~ 2.5crDD probing the cavity size (Fig. 
CH middle). 

Dynamic arrest of big particles is essentially driven 
by soft-sphere repulsion, as it would be for a sys- 
tem without size disparity, and a consistent descrip- 
tion of dynamic features is achieved within the stan- 
dard A2-scenario of the MCT. Still, it must be men- 
tioned that the obtained value of the exponent pa- 
rameter A = 0.887 is signifficantly higher than typ- 
ical values A ^ 0.7 usually obtained within the A2- 
scenario |^ y, |j, p . This high value of A can be in- 
terpreted as a signature of the depletion mechanism 
induced by neighboring small particles, which leads 
to the observed weak clusterization of big particles. 
Though only weakly competing with soft-sphere re- 
pulsion, the depletion mechanism yields precursor ef- 
fects of an incoming higher-order scenario (A — > 1), as 
strong stretching for the long-time decay in _F'aa('?, i) 
(Fig [3d) , or sublinear behavior in the mean squared 
displacement (Fig Ob). These effects are weaker for 




FIG. 11: (color online) Typical slabs (of thickness 5(tdd) 
of different configurations at T = 0.30. From top to down 
a:smaii = 0.1, 0.6, and 0.8. Dark: A- and B-particIes. Light: 
C- and D-particles. 



the composition a;smaii — 0.3, yielding a lower expo- 
nent parameter A = 0.880. 

b) High concentration of small particles 

For composition a;smaii = 0.8, the population of 
small particles is sufficient to provide an efficient de- 
pletion mechanism, yielding a strong clusterization 
of the big particles (Fig. ^2 bottom). Competi- 
tion between soft-sphere repulsion and depletion leads 
to anomalous relaxation features for the big parti- 
cles. Hence, Faa('?,0 shows logarithmic decay for 
q ~ 5.6CTQQ (Fig. |S^), i.e., for wavelengths ^ ctdd 
probing the size of the small particles. 
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FIG. 12: (color online) Radial distribution function for 
D-C and D-D pairs, for Xamaii = 0.1 and T = 0.50. 



Though depletion effects are evidenced by cluster- 
ization of big particles, it is worth emphasizing that 
a direct comparison with results for the effective one- 
component systems dis playing the higher-order MCT 
scenario |al3)lEliill3jll3jllJ| cannot be made. In the 
effective one-component systems the small particles 
are absent and the depletion mechanism is described 
by an effective short-ranged attraction between the 
big particles [l9j. However, concerning dynamics, the 
validity of this approximation does not only require a 
much smaller size, but also a much smaller mass for 
the small particles [38|. In the case here investigated 
big and small particles have the same mass. There- 
fore, a proper proof of the existence of an underly- 
ing higher-order MCT scenario can only be provided 
by solving the MCT equations for the two-component 
system. 

Density-density correlators for small particles at 
a^smaii — 0.8 exhibit a complex decay. An apparent 
logarithmic decay occurs at long times (Fig. ^). 
For T = 0.30 the onset of this decay occurs at 
t K, 60, which corresponds to mean displacements 
((ArD)^)^/^ ^ 2crDD (Fig. Et). If logarithmic decay 
is again interpreted as the result of a competition be- 
tween bulk-like dynamics and confinement, the former 
distance would correspond to the length scale beyond 
which confinement induced by the big particles affects 
dynamic arrest of the small particles. For smaller 
length scales, the only relevant arrest mechanism is 
bulk-like caging induced by neighboring small parti- 
cles. Hence, for i ^ 60 small particles do not exhibit 
clear signatures of anomalous relaxation, as suggested 
by Figs. Et andEt. 

c) Low concentration of small particles 

For Xsmaii = 0.1, most of the small particles are 
isolated in the slow matrix of big particles (Fig. 1111 
top). However, there is also a signifficant population 
of pairs of neighboring small particles. We have esti- 
mated that about a 7 % of the small particles have a 
neighboring small particle within a distance I.Sctdd- 
Clusters of three or more neighboring small particles 
for this inter-particle distance are very rare. 

The existence of pairs of neighboring small particles 
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FIG. 13: (color online) Symbols: For composition a^snmii = 
0.1, mean squared displacement (panel (a)), and density- 
density correlator (panel (b)) for the subset of D-particles 
which initially have at least one neighboring small particle 
within a distance r < I.Sctdd. Panel (a) shows data for 
T — 0.30 and 0.40. Dashed curves are the corresponding 
data by averaging over all the D-particles. Solid lines in- 
dicate linear or sublinear behavior. Panel (b) shows data 
for different g-values at T = 0.30. Dashed curves are data 
by averaging over all the D-particles (from top to down 
go-DD = 0.92, 2.00, 3.23, and 5.02). The solid line indi- 
cates logarithmic relaxation. 



yields a sharp peak at distances ctdd S f ^ S.Sctdd for 
the radial distribution function for D-C (5dc('')) and 
D-D {griD{r)) pairs, as shown in Fig. ^1 Hence, even 
for so low populations of small particles, it may be 
expected that competition between bulk-like dynam- 
ics and confinement occurs for a signifficant fraction 
of particles. In order to test this hypothesis we have 
calculated mean squared displacements and density- 
density correlators for the subset of D-particles which 
initially (i.e., at the time origin for the calculation of 
both quantities) have at least one neighboring small 
(C or D) particle within a distance r < I.Sctdd- This 
distance approximately corresponds to the location of 
the maximum in 5dc(^) and guuif)- Fig. 1131 shows 
((Arn)^) and Fdd((7, t), computed for the former sub- 
set of D-particles, and compared with the correspond- 
ing quantities computed for all the D-particles (pre- 
viously shown in Figs. ^ and 1 10b). Relaxation of 
D-particles having neighboring small particles is, at 
intermediate times, faster than the average over all 
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D-particles, which is approached only for very long 
times corresponding to the onset of the linear diffu- 
sive regime. 

Differently from the average over all D-particles, 
the subset initially having neighboring small parti- 
cles does exhibit anomalous relaxation features. Mean 
squared displacements display sublinear behavior over 
two time decades. The corresponding exponent de- 
creases with decreasing temperature (Fig. \Y6h). 
Density-density correlators show a concave-to-convex 
crossover (Fig. 113b). Logarithmic relaxation is ob- 
served for wavevectors q ~ 2.5ctqq (Fig. Et), i-c, for 
length scales ~ 2.5(Tdd which probe the first minimum 
in gucir) and guuir). 

We have also computed ((ArD)^} and Fdd(<7, i) 
for the subset of D-particlcs which initially do not 
have any neighboring small particle within a distance 
r < 4.5crDD, i-e., for initially isolated D-particles. Re- 
sults for this subset (not shown) are hardly distin- 
guishable from those averaged over all D-particles and 
reported in Figs. 0^ and llOb . Hence, initially iso- 
lated D-particles do not show anomalous relaxation 
features as sublinear behavior for mean squared dis- 
placements or concave-to-convex and logarithmic de- 
cay for density-density correlators. 

d) Very large size disparity 
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In this subsection we show and discuss results for 
the small particles, for a very large size disparity {S = 
8.0) at a single composition a;smaii — 0.6. Hence, we 
select a composition for which small particles exhibit 
anomalous relaxation features at disparity S = 2.5 
(Figs. ^,03 andinj. We also use the same packing 
fraction, (p — 0.53, that for the case S — 2.5. We ob- 
serve a relaxation scenario rather different from that 
obtained for disparity 6 = 2.5. Fig. 114b displays 
density-density correlators for A-A and D-D pairs, at 
a wavelength of about 10 times the size of the small 
particles. In order to provide a clearer visualization of 
the decay for the latters, data are plotted as a func- 
tion of tT^/^, i.e., rescaling the time by the thermal 
velocity. As expected, big particles exhibit a much 
slower dynamics. Correlators for small particles do 
not show a slow decay at any temperature. Decorre- 
lation occurs in an essentially exponential way down 
to very small values of FjjD{q, t), where a tail of small 
amplitude arises and extends to very long times. 

Below some given temperature a background clearly 
distinguishable from zero arises. The amplitude of the 
background increases as the system is cooled down. 
Qualitatively similar results (not shown) are observed 
for other wavevectors. We identify the background 
amplitude as an operational non-ergodicity parameter 
that increases from some critical value f^. By visual 
inspection of Fig. I14h . it is difficult to establish wether 
fg defined in this way is zero, but it is evident from 
data at T = 0.40 that it takes, as much, a extremely 
small value. This behavior resembles features charac- 
teristic of MCT transitions of the type-A, which are 
defined by a zero value of the critical non-ergodicity 
parameter, in contrast to the finite value defining stan- 
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FIG. 14: (color online) Simulation results at composition 
a^smaii = 0.6 and size disparity 5 = 8.0. Symbols in pan- 
els (a) and (b) are data for, respectively, Foo{q,t) and 
Fli{q,t). Symbols in panel (c) are data for the mean 
squared displacement of D-particles divided by the time 
t. Thick solid (T = 0.03) and dashed (T = 0.20) curves 
in the three panels are the corresponding results for A- 
particles. Horizontal straight lines in panels (a) and (b) 
indicate the zero value. 



dard type-B transitions. A recent realization of this 
scenario has been reported for a system of dumbbell 
molecules [39|, |40|, |41| . While for moderate molecular 
elongations a standard relaxation scenario is observed, 
for small elongations angular correlators exhibit fea- 
tures analogous to those of Fig. 114b l4ll. Theoreti- 
cal calculations for this latter system |39| relate such 
features to the existence of a MCT transition of the 
type-A. 
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Fig. inb displays self-correlators, F^{q,t) = 
Ej expjiq • [ra,j{t) — ra,j{0)]}, for A- and D-particles. 
As observed for Faa(9,0; self-correlations for A- 
particles show a much slower dynamics. Self- 
correlators for D-particles display a striking result. 
As for density-density correlations (Fig. [Tlk). they 
exhibit a fast decay followed by a long tail. How- 
ever, they decay to zero at all the investigated tem- 
peratures, even at those where -Fdd('?, t) show a finite 
long-time limit. Fig I14h shows results for the mean 
squared displacement. Data are divided by t in order 
to evidence that D-particles reach the diffusive regime 
(horizontal limit in this representation) at long times. 
Data in Figs. 114b and 1 14b show that the self-motion of 
small particles is ergodic at all the investigated tem- 
peratures. A small particle can reach regions arbi- 
trarely far from its initial position. However, accord- 
ing to results in Fig. 114b . coherent dynamics are non 
ergodic below some given temperature. Hence, two 
small particles can explore the whole structure of the 
matrix of big particles, but they will never fully decor- 
relate from each other, i.e., the distance between them 
will always be finite. 

These features are consistent with early MCT theo- 
retical calculations by Bosse and Thakur |42l |43| for a 
binary mixture of hard spheres. Though in that work 
no information is given about relaxation features, cal- 
culations are reported for non-ergodicity parameters 
of self- (fg) and density-density (/,) correlators for 
the small particles. For a symmetric c omp osition and 
large disparity size (S ~ 10) it is found |43| that /| = 
at packing fractions where fq > 0, in agreement with 
results reported here with temperature as the control 
parameter. 

e) A possible unified picture 

The whole picture here reported for the small par- 
ticles at large {S — 2.5) and very large {S = 8.0) size 
disparity provides a connection with MCT theoretical 
results by Krakoviack in Refs. |2a, |23 for a mixture 
of fixed and mobile hard spheres. In that work the 
dynamic phase diagram displays an A- and a B-line 
in the plane 0fix-0mob, where (pax and 4>moh are the 
packing fractions of, respectively, the fixed and mobile 
particles. The A- and B-lines merge at a higher-order 
point (namely an As-point). The B-line extends from 
the As-point to the limit (pRx — 0, where the liquid 
of hard spheres is recovered. The A-line extends from 
the Aa-point to the Lorentz gas limit at ipmoh — 0. 

For high concentrations of the mobile particles, the 
matrix of fixed particles does not yield signifRcant con- 
finement effects. The only caging mechanism is nor- 
mal hard-sphere repulsion at short length and time 
scales, and dynamic correlators exhibit a standard 
two-step decay [23 • The transition point is of the 
type-B and hence the jump of the lon g-ti me limit of 
the density-density correlator is finite [23, providing 
a non-zero value of the critical non-ergodicity param- 
eter. 

For high dilution of the mobile particles, hard- 
sphere repulsion at short length and time scales can- 



not yield temporary caging, and density-density cor- 
relators display a fast decay to values close to zero 
J23|. At longer time scales the mobile particles probe 
the structure of the confining matrix of fixed parti- 
cles, which leads to a "mesoscopic" caging, character- 
ized by a length scale larger than that characteristic 
of bulk-like hard-sphere repulsion. As a consequence 
of this large-scale caging mechanism (confinement), 
density-density correlators for the mobile particles ex- 
hibit a long tail of small amplitude after the fast mi- 
croscopic decay |27|. At the transition point, of the 
type-A, the long-time limit does not exhibit a finite 
jump but grows up continuously !27|. providing a zero 
value for the critical non-ergodicity parameter. 

As mentioned above, at moderate concentrations of 
fixed and mobile particles, a higher-order Aa-point 
arises (2^] as a consequence of the competition be- 
tween the mentioned "microscopic" and "mesoscopic" 
caging mechanisms. Relaxation features have not 
been reported in Refs. |23,|27J for state poin ts close to 
the As-point. However, as stressed in [23, they will 
necessarily display the anomalies reported here, as a 
mathematical consequence of the value A = 1 defining 
the Aa-point. 

How do results in Refs. |2(1 l27l | compare with re- 
laxation features presented here?. First it is worth 
emphasizing that, differently from Krakoviack's work, 
high dilution of the small particles is not the key in- 
gredient for yielding a type-A transition for the latters 
in the present system. Indeed, data reported here for 
S = 2.5 and a;sman = 0.1 show features rather differ- 
ent from those characterizing type-A transitions. It 
must be noted that, for these control parameters, the 
system here investigated is much denser than at state 
points close to the A-line for the mixture of fixed and 
mobile particles |2a|. Despite its low density, diffu- 
sion of the mobile particles in the latter system is 
blocked at large length scales due to the absence of 
percolating free volume. In the system here investi- 
gated, small particles can diffuse at higher densities 
due to the non-static nature of the confining matrix. 
The slow motion of the big particles creates regions of 
sufficient local free volume which facilitate diffusion 
of the small particles. As a consequence of high den- 
sity, short-range bulk-like caging is a relevant arrest 
mechanism for small particles and leads to a slow de- 
cay of dynamic correlators (Fig. llOb^. different from 
features characteristic of type-A transitions. Hence, 
for disparity S = 2.5, type-A transitions cannot exist 
at any mixture composition. 

A way for yielding a type-A scenario for the small 
particles in the present system, at a fixed packing frac- 
tion and mixture composition, is by increasing the 
size disparity, as shown above. In this situation small 
particles move in a medium of low local density, and 
short-range bulk-like caging is supressed. This effect 
is illustrated in Fig. E| by displaying a typical slab 
of the simulation box. However, small particles are 
blocked at large length scales below some given tem- 
perature where the tail of the soft-sphere potential is 
probed. At that temperature, despite of local vibra- 
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FIG. 15: (color online) Typical slab (of thickness IOctdd) 
a,t T — 0.03 and aJsmaii ~ 0.6 for size disparity 5 = 8.0. 
Dark: A- and B-particles. Light: C- and D-particles. 



tions of the confining matrix, the effective radius of 
the big particles becomes so large that the percolat- 
ing free volume vanishes, leading to a dynamic picture 
analogous to the type-A scenario observed for the sys- 
tem of Refs. miia. 

Relaxation features observed for the small particles 
at S = 8.0 and Xsmaii = 0.6 suggest the presence of 
a nearby A-line, originating from large-scale caging 
(confinement) induced by the slow matrix of big par- 
ticles. For small size disparity, 6—^1, one obviously 
recovers the standard MCT scenario predicted for the 
monoatomic liquid of har d sp heres and widely ob- 
served in mixtures of hard [25| , soft |4J| , or Lennard- 
Jones spheres H, where small disparity is introduced 
just to avoid crystallization. The standard scenario, 
originating from short-range bulk-like caging, is char- 
acterized by a nearby B-line. 

Increasing the size disparity at fixed composition 
a^smaii = 0.6 will weaken the effects of bulk-like caging 
and strengthen those associated to confinement. This 
feature suggests that a crossover from a B- to an A-line 
will occur by increasing the value of (5, in analogy with 
results for the system investigated in Refs. |23, 123 , 
where such a crossover is obtained by varying the mix- 
ture composition. For the latter system, the B- and A- 
lines merge at a higher-order point. The resemblance 
of results reported here for S = 2.5 and Xsmaii = 0.6 
(Figs. 0Jdi Eb and EJ with relaxation features char- 
acterizing higher-order MCT transitions suggest an 
analogous merging. The former results would origi- 
nate from the presence of a nearby B-line (due to the 
finite value of the non-ergodicity parameters), end- 
ing at a nearby higher-order point (which would pro- 
duce the observed anomalous relaxation features), or 
at least at an A2-point with A very close to unity. 

VI. CONCLUSIONS 

We have carried out simulations on a mixture of big 
and small particles. Slow dynamics have been inves- 



tigated for a broad range of temperature and mixture 
composition. The introduction of a signifBcant size 
disparity yields very different time scales for big and 
small particles, reaching differences of 2-3 decades in 
diffusivity for the lowest investigated temperatures. 
This model exhibits non-conventional relaxation fea- 
tures. Mean squared displacements display sublin- 
ear behavior at intermediate times. The exponent 
for the corresponding power law decreases by decreas- 
ing temperature. By varying temperature or wavevec- 
tor, a concave-to-convex crossover is obtained for the 
shape of the decay of density-density correlators. At 
some intermediate point of this crossover, the decay 
is purely logarithmic. 

These anomalous relaxation features, which are 
observed over time intervals extending up to four 
decades, strongly resemble predictions of the Mode 
Coupling Theory (MCT) for state points close to 
higher-order transitions, which originate from the 
competition between different mechanisms for dy- 
namic arrest. By varying the mixture composition, 
anomalous relaxation is displayed by both the big 
and the small particles. For the big particles we sug- 
gest competition between soft-sphere repulsion and 
depletion effects induced by neighboring small par- 
ticles. For the small particles we suggest competition 
between bulk-like dynamics and confinement, respec- 
tively induced by neighboring small particles and by 
the slow matrix of big particles. 

We have also performed simulations at a fixed com- 
position for a very large size disparity. A new relax- 
ation scenario arises for the small particles, showing 
features characteristic of nearby MCT transitions of 
the type-A. This feature provides a connection with 
MCT theoretical results for a mixture of mobile and 
static particles, which report a dynamic phase dia- 
gram displaying an A- and a B-line merging at a 
higher-order point. A similar crossover is suggested 
for the system here investigated by varying the size 
disparity. 

Simulation results reported here do not constitute a 
rigorous proof of an inherent MCT scenario as the one 
described above. A proper answer to this question can 
only be provided by solving the corresponding MCT 
equations. However, the highly non-trivial observed 
analogies suggest to consider it as a plausible hypoth- 
esis. Though being beyond the scope of this article, 
the observation of confinement effects also suggests to 
test a more speculative hypothesis: a MCT scenario 
for a liquid of spheres in fractal dimensions as the ori- 
gin of the observed relaxation features for the small 
particles. To our knowledge no MCT calculations of 
mean squared displacements or dynamic correlators 
are available for the latter or similar systems. 
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